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ABSTRACT 



The magnetic field in the local interstellar medium (ISM) provides a key indi- 
cator of the galactic environment of the Sun and influences the shape of the helio- 
sphere. We have studied the interstellar magnetic field (ISMF) in the solar vicinity 
using polarized starlight for stars within 40 parsecs of the Sun and 90° of the helio- 
sphere nose. In Frisch et al. (2010, Paper I) we developed a method for determining 
the local ISMF direction by finding the best match to a group of interstellar polariza- 
tion position angles obtained towards nearby stars, based on the assumption that the 
polarization is parallel to the ISMF In this paper we extend the analysis by utilizing 
weighted fits to the position angles and by including new observations acquired for 
this study. We find that the local ISMF is pointed towards the galactic coordinates 
£,b= 47° ± 20°, 25° ± 20°. This direction is close to the direction of the ISMF that 
shapes the heliosphere, £,b— 33° ±4°, 55° ±4°, as traced by the center of the "Ribbon" 
of energetic neutral atoms discovered by the Interstellar Boundary Explorer (IBEX) 
mission. Both the magnetic field direction and the kinematics of the local ISM are 
consistent with a scenario where the local ISM is a fragment of the Loop I superbub- 
ble. A nearby ordered component of the local ISMF has been identified, in the region 
£m 0° — > 80° and bm 0° — > 30°, where PlanetPol data show a distance-dependent 
increase of polarization strength. The ordered component extends to within 8 parsecs 
of the Sun and implies a weak curvature in the nearby ISMF of ~ 0.25° per parsec. 
This conclusion is conditioned on the small sample of stars available for defining this 
rotation. Variations from the ordered component suggest a turbulent component of 
~ 23°. The ordered component and standard relations between polarization, color ex- 
cess, and H° column density predict a reasonable increase of iV(H) with distance in the 
local ISM. The similarity of the ISMF directions traced by the polarizations, the IBEX 
Ribbon, and pulsars inside the Local Bubble in the third galactic quadrant suggest that 
the ISMF is relatively uniform over spatial scales of 8-200 parsecs and is more similar 
to interarm than spiral-arm magnetic fields. The ISMF direction from the polariza- 
tion data is also consistent with small-scale spatial asymmetries detected in GeV-TeV 
cosmic rays with a galactic origin. The peculiar geometrical relation found earlier be- 
tween the cosmic microwave background dipole moment, the heliosphere nose, and 
the ISMF direction, is supported by this study. The interstellar radiation field at ~ 975 
A does not appear to play a role in grain alignment for the low density ISM studied 
here. 

Subject headings: ISM: magnetic field, clouds, solar system: general, Galaxy: solar 
neighborhood, cosmology: cosmic microwave background 
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1. Introduction 



The characteristics of the interstellar magnetic field (ISMF) in the warm, low density, partially 
ionized interstellar gas near the Sun are difficult to study. A new diagnostic of the magnetic field 
in the Local Interstellar Cloud (LIC) surrounding the helio sphere is provided by the giant arc, 
or 'Ribbon', of ener getic neutral atoms (ENAs) discovered by the Interstellar Boundary Explorer 
(IBEX) spacecraft (IMcComas et al.ll2009l . l201ll) . ENAs are formed by charge-exchange between 
interstellar neutral hydrogen atoms and the solar wind and other heliospheric ions. The ENA 
Ribbon is observed towards sightlines that are perpendicular to the ISMF direction as it drapes over 
the heliosphere, and is superimposed on a distributed ENA flux that traces the global heliosphere 



(App endix IA1 provides additional information on the Ribbon and Ribbon models). In Frisch et al 



(|2010al Paper I) we showed that the measurements of interstellar polarization towards nearby stars 
yields a best-fit ISMF direction that is close to the field direction obtained from the IBEX Ribbon. 
In this paper, we further study the connection between the magnetic field sampled by the IBEX 
Ribbon and the magnetic field in the local ISM. 

Starlight polarized by magnetically aligned dust grains in the ISM shows that the magnetic 
field in our galactic neighborhood has two large-scale components, a uniform field parallel to the 
galactic plane that is directed towards £~ 82. 8111 and extends beyond a kiloparsec, and a giant mag- 
netic structure, Loop I, that is within 100 pc and dominates the northern sk y and probably reaches to 
the solar location. The ISMF of Loop I is traced by opt ical polarizations ( Mathewson & Fordll970[ 
Heil esl 19761: Santos et aDboi ll i feerdyugin et alJboi l|) . t he polarized radio continuum (iBerkhuijsenl 



1973 



Wolleben 



2007|) . and Faraday rotation measures (|Taylor et al.ll2009b iMao et al. 



201dh . It is 



a superbubble fo rmed by stellar winds and supernova in the Sco-Cen A ssociation during the past 
~ 15 Myrs (e.g. Ide Geus1ll992l : iFrischll 19951 fl99rl lHeileslll998aL I2009L The interstellar magnetic 
field is swept up during superbubble e xpansion, creating a magnetic bubble that persists through 
the late stages of the bubble evolution (|Tilley et al.ll2006l) . Optical polarization and reddening data 
show that the eastern parts of Loop I, 1= 3° — 60°, b> 0°, are within 60-80 parsecs of the Sun 
(ISantos et al.ll201ll : LFrisch et al.ll201l[) . If Loop I is a spherical fea t ure, the ISMF near the Sun 
should be associated with the Loop I superbubble (e.g. Frisch 1990 : Heiles 1998b a|). The mag- 
netic structure of Loop I provides a framework for understanding the magnetic field in the local 
ISM and the relation between the galactic magnetic field and the magnetic field that shapes the 
heliosphere. 

In this paper we use high- sensitivity observations of the polarization toward nearby stars as 
a tool for determining the direction of the ISMF near the Sun, within ~ 40 parsecs. This study 



'in this paper, galactic coordinates are denoted by £, b, and ecliptic coordinates by A, j3. Plots in ecliptic coordi- 
nates are labeled as such. 
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builds on our earlier study (Paper I) which presents a new method for determining the ISMF di- 
rection, based on testing for the best match between polarization position angles and the ISMF, 
where we found that the best-fit local ISMF direction in the galactic hemisphere regions agrees, to 
within the uncertainties, with the ISMF direction defined by the center of the IBEX Ribbon arc. 
Our new study differs from the first study in that the method is broadened to include weighted fits 
(©, which is possible because of the inclusion of new data (©. The best-fit ISMF direction found 
from both unweighted fits ( ^4.11 ) and weighted fits ( 514.21 ) is close to th e direction of the ISMF 
shaping the heliosphere ( ^5.1| ). We also look closely at the PlanetPol data ( Bailey et alJboiO , S15.3I ) 
and use the upper envelope of the polarization-distance curve ( §5.3.11 Appendix to show that 
there is an ordered component in the local ISMF with a direction that varies weakly with distance 
and extends to within 8 parsecs of the Sun ( §5.3.21) . The ordered component provides information 
on magnetic turbulence. An anti-correlation of polarization with radiation fluxes suggests ineffec- 
tive radiative torques ( ^5.3.3[ Appendix |D1). The best-fit ISMF direction is compared with the field 
direction derived from the IBEX Ribbon, and the relation between the LIC and G-clouds ( §5.11) . 
For our new ISMF direction and the new IBEX heliosphere nose direction, the symmetry of the 
cosmic microwave background (CMB) dipole moment with respect to the heliosphere nose direc- 
tion is still apparent ( ^5.3.4| ). The similar directions of the ISMF over scales of several hundred 
parsecs suggests an interarm-type field ( §5.4.11) . The local ISMF direction also appears to affect 
anisotropics observed in GeV-TeV galactic cosmic rays ( §5.4.21) . Results are summarized in §|6] 



2. Method of finding the best-fit ISMF direction 



Optical polarization is a n important probe of the di rection of the ISMF in in terstellar clouds , 
both oyer large spatial scales (iMathewson & Fordlll970|) and the very local ISM (!Tinbergenlll982 : 
Frischlll990[ iFrisch et alJboiOah . The observed pola rization from aligned grains depends on the 
total column of dust and the grain size distribution (|Mathislll986|) . the fraction of asymmetrical 
dust grains and alignment efficiency of the grains (Whittet et al. 2008), and the projection of 
the magn etic field onto the sky. The axis of lowe s t average grain opacity is para llel to the ISMF 
direction ( Davis & Greenstein 1951 : Martin 1971 : Roberge 2004 : Lazarian 2007|) . @ Polarization 
is a pseudo-vector entity (always in the range 0° — 180°), and in denser regions radiative transfer 
effects can influence the observed polarization both due to turbulence in the probed medium (Jones, 
Dickey & Klebe 1992; Ostriker, Stone & Gammie, 2001), and due to the influence of foreground 



2 Leverett Davis provided the first viable model for the relation between the ISMF and starlight polarization, and 
also first proposed the existence of the solar wind cavity in the ISM, now known as the heliosphere. He proposed a 
solar wind ca vity of radiu s ~ 200 AU partly out of the necessity of excluding the interstellar magnetic field from the 
solar system ( Davisfl955 ). 



-5- 



depolarizing screens (|Andersson & Potter! l2006t I Anders sonl |2 1 2|) . ISM within 10 parsecs has 
low average ISM densities, iV(H)< 10 18 cm -2 and < n >~ 0.1 cm~ 3 , suggesting there is very 
little nearby dust. Therefore depolarizing foreground screens and collisional disruption of grain 
alignment are les s likely locally com pared to distant dense cloud regions. This low density ISM is 
partially ionized (|Frisch et al.ll201l|) so that gas and the ISMF will be tightly coupled. 



Two properties of polarized starlight could in principle be exploited to obtain the ISMF di- 
rection causing the polarization: the position angle of the polarizations, which are measured with 
respect to a great circle meridian and increase positively to the east, and the polarization strengths 
that reach a maximum for sightlines perpendicular to the ISMF for a uniform medium. The draw- 
backs with using polarization strengths to determine the ISMF direction are that polarization varies 
with the dust column density of the sightline, and polarizations are weaker near the magnetic pole 
where the statistical significance of detections will be lower. Polarization directions are insensitive 
to the dust column density, but sensitive to magnetic turbulence that could distort the field over 
either small or large scales. 

Paper I introduces a new method for determining the ISMF direction based on finding the 
best fit to an ensemble of polarization position angles. The analysis makes the assumptio n that the 
polari zation position angle is parallel to the direction of the interstellar magnetic field (Roberge 
2004|) . This assumption is justified for the diffuse clouds in Loop I, where the field directions 
determined from optical polarization data and synchrotron emission in high-latitude regions agree 
(ISpoelstralll972|) . The ensemble of data used in the analysis are chosen to meet a set of selection 
criteria. These data and the criteria are discussed in the following section (©. All measurements 
meeting the criteria are used in the analysis, including multiple measurements of the same star by 
different observers. 

The ISMF direction is derived by assuming that a single large-scale ISMF close to the Sun 
is aligning the grains, so that the polarization direction would then be parallel to a great-circle 
meridian of the "true" magnetic field direction. The best-fit ISMF is calculated to be the direction 
that corresponds to the minimum value of: 



Fi = F(Bi) 



N" 



n=l 



sin(g n (^i)) 
G n 



(1) 



where 9 n (B{) is the polarization position angle PA Q for star n, calculated with respect to the i th 
possible interstellar magnetic field direction £>;, and where the sum is over N stars. G n is the 
weighting factor for each star n based on the measurement uncertainties of the polarization position 
angle for that star. A grid resolution of 1° in longitude and latitude is used for the spacing of 
Bi values. An alternate approach to find the best-fit field direction -Bi=b es t that maximizes the 
cosecants, i.e. maximize \G n / sin(9 n (B{)\ , generated numerical glitches related to dividing by 
zero, and is not used. 
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Paper I calculated the ISMF direction with the weighting factor G n =l . Statistical weights were 
ignored since the northern hemisphere data included a large number of recent high-sensitivity data 
collected at ppm accuracy, while the southern hemisphere data were dominated by data collected 
in the 1970's that are less sensitive (la >60 ppm, Weighting data points by the uncertainties 
yielded a magnetic field direction biased by the difference in uncertainties in the two hemispheres. 
The best-fit ISMF direction from Paper I is listed in Tabled] In this paper, new data are added 
and the unweighted fit is recalculated ( §4.1| ). 



New high-precision polarization data make it possible to evaluate the best-fit ISMF direction 
by utilizing all of the position angles in the sample, but weighting individual values according to 
their statistical significance. By combining data based on weighted values, the emphasis is implic- 
itly placed on the ISMF field in the distant portions of the sampled ISM, where polarizations may 
be stronger relative to the measurement errors (depending on the patchiness of the ISM, §5.31 ). The 
probability distribution for position angles does not have a normal distribution for weak polariza 



tions, P /dP< 6. The weighting function is based on the results of iNaghizadeh-Khouei & Clarke 



(I1993a|) that give the probability distribution for the position angle of linear polarization. For true 
polarization p Q with a measurement uncertainty of a, and the true position angle in equatorial co- 
ordinates 8 , the probability G r n (6' b s ; 9 Q , P Q ) of observing a position angle o b s (which will be 0° 
for the "true" value of a perfectly aligned PA) is: 

G n (9 ohs ; o , P a ) = -L {4= + exp(r /o 2 ) [1 + erf (77 )]}exp(-^) (2) 



where r] = ^ cos [2(8 ohs — o )) and P Q = ^f, and the function erf(z) is the Gaussian error 

function erf(Z) = ^= J Q Z exp(—t 2 ) dt. The goal is then to choose a procedure that maximizes the 
contributions of high-significance (e.g. high G n ) measurements, and minimizes the contributions 
from data points with large values of sin(0 n j). When -Bi=b es t is calculated using weighted fits 
( §4.21) . the function G n in eqn. |2]is adopted as the weighting factor for each star n in eqn. [Q 



3. Polarization data used in study 

The Local Bubble void appears as an absence of both interstellar gas and dust within ~ 70 
parsecs, particularly in the third and fourth galactic quadrants. Figure \T\ shows the distribution of 
interstellar dust within 100 pc, as traced by color excess _E(B-V)o The interior of the Loop I su- 
perbubble appears as a void in the distribution of interstellar dust within 100 parsecs in the fourth 



Stars with SV > 0.06 mag, according to the Hipparcos variability index H6, are omitted. 
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galactic quadrant (£> 270°). Models of Loop I as a spherical object place the Sun in or close to 
the rim of the shell ( ^5.21) . This combination of properties is the basis for selecting polarized stars 
for this analysis within 90° of the heliosphere nose and 40 parsecs of the Sun. The heliosphere 
nose is remarkably close to the direction of the galactic center. IBEX-LO measurements of the 
23.2 ± 0.3 km s -1 flow of interstell ar He° through t he heliosphere provide the nose direction of 
£= 5.25°±0.24°,6= 12.03°±0.51° dlVIcComas et alJ2012h . Most interstellar gas within 40 parsecs 
is also clumped within 20 parsecs of the Sun, as shown by interstellar H° and D° (Figure |2]> and 
the fact that the average number of interstellar Fe + , Ca + , and Mg + velocity components is distant- 
indepen dent for stars withi n 30 parsecs, and rises only slightly (~ 20%) for stars that are 30-40 
parsecs (|Frisch et al.ll201l|) . This spatial interval provides for a well-defined region for study, ac 



cepts most stars in the region of Tinbergen's origina l "patch" of nearby polarization (|Tinbergen 



19821 Paper I), includes the region where PlanetPol (|Bailey et al.ll2010l) obtained v ery high sen 



sitivity polarization data ($3, and includes nearby stars observed towards Loop I by iSantos et al 
J201lh . 



These data show that the polarizations in the southern galactic hemisphere tend to be larger 
than polarizations in the northern galactic hemisphere. For the special case where the star sample 
is restricted to stars meeting the spatial and spectral criteria above, and to recent higher precision 
data, then the polarizations of southern stars with P/dP> 2.5, are 0.024% ± 0.017%, while for 
northern stars they are 0.003% ± 0.005%. Whether this difference, originally noted by Tinbergen 
(1982), is accurate will require additional high-sensitivity data (where P/dP gives the polarization 
strength divided by the la uncertainty of the polarization). 

New polarization measurements were also collected for this study using telescopes of the 
Laboratorio Nacional de Astropfisica (LNA) in Brazil and the Nordic Optical telescope (NOT) 
on the islan d of La Palma in the C anary Islands. The LNA data were acquired with the CCD 
polarimeter (Magalh aes et al.lll996l) on the 0.6m telescope at Pico de Dios, during the months of 
June 2008, August 2010, September 2010, and May 2 011. The NOT data were acquired using the 
TurPol polarimeter (|Piirolal 1 1 973c butters et al.ll2009|) on the 2.52m telescope on the island of La 
Palma in the Canary Islands during June 2010. Optical polarimetry is a differential measurement 
and good observing conditions are not as critical as for high accuracy photometry. ISM polarization 
data also sample a relatively constant polarization source, allowing repeated observations over 
different days or from different sites to maximize sensitivity and minimize systematic uncertainties. 
These new data are listed in Table |2] The columns in this table are HD number, galactic longitude 
and latitude, star distance, the polarization in units of 10~ 5 , and the polarization position in the 
equatorial coordinate system. The last column gives the observatory where the data were acquired. 



The CCD imaging polarimeter at Pico dos Dias provided polarizations in the B-band. The 
instrument acquires CCD images of the orthogonal polarizations, in sets of 8 waveplate positions, 



-8- 



with several images at each plate position. In the case of constant polarizations such as interstellar 
polarizations, individual observations can be combined to reduce uncertainties. Stars brighter than 
~ 7 mag were observed through a neutral density filter. 

The TurPol data were obtained using the UBV mode of the polarimeter with high voltage 
for the R and I photomultipliers switched off. The effective wavelengths of the UBV bands are 
3600, 4400, and 5300 A, respectively. With this mode, typical count rates of 10 6 counts s _1 were 
achieved giving a precision of 0.01% for a 3a detection. Observations of two zero-polarization 
standards were acquired each night at the same precision in order to determine the instrumental 
polarization. The position angle calibration was obtained with observations of two high polariza- 
tion standard stars. For interstellar polarizations, position angles do not depend significantly on 
the wavelength, allowing the use of weighted-averages from the UBV bands to achieve the highest 
possible precision; these 'broad-band' position angle measurements are listed in Table |2] 



The early iTinbergenl (|1982|) survey, made during 1973-1974 from La Silla (Chile) and Har- 
tebeespoortdam (S. Africa) in the southern hemisphere, and Leiden (Netherlands) in the northern 
hemisphere, remains the most extensive all-sky survey of weak polarizations, with la sensitivity 
on the Q and U Stokes parameters of 0.006-0.009 %, depending on the instrument and data set. 
These weak polari zations are biased towards larger pola rization strengths, since P 2 = Q 2 + U 2 
is always positive. iNaghizadeh-Khouei & Clarke! (Il993bh evaluated the statistics of stars observed 
from each observing station in the Tinbergen data set, and concluded that residual instrumental po- 
larizations may remain in the K-star data collected from Hartebeespoortdam. We have examined 
the data from Hartebeespoortdam for non-active K-stars that also have P/dP> 2.5. The seven 
stars that match the data criteria are located at angles of 86° ± 32° from the ISMF direction defined 
by the center of the IBEX Ribbon arc. If the direction of the ISMF at the heliosphere is the same 
as in the local ISMF direction, then for a uniform ISM the maximum polarizations are expected at 
angles of 90° from the ISMF pole. We therefore assume that any contribution of instrumental po- 
larization to the Hartebeespoortdam data is insignificant, while acknowledging that high-sensitivity 
observations and temporal monitoring of the polarizations of these stars are needed. 



This study also includes the IBailey et alj (|2010l) polarizations of nearby stars obtained with 
the PlanetPol instrument at the William Herschel telescope. These linear polarization data were 
collected over a broad red band (maximum sensitivities at 7000-8000A), and achieved typical 
accuracies of parts-per-million. PlanetPol data were also included in the Paper I analysis, but 
without screening for stellar activity. Twenty-six stars in the PlanetPol catalog that passed the 
selection criteria abov e were incorporate d into our analysis of the local ISMF direction and are 
listed in Appendix IB] IBailey et al.l (|2010 ) measured 18 northern hemisphere stars that were also 
observed by Tinbergen, and detected weak polarizations towards roughly half of the stars that were 
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unpolarized at the sensitivities of the Tinbergen data0 PlanetPol data showed that polarizations in 
the RA > 17 H region increased with distance, and measured a weak interstellar polarization of 
17.2 ± 1.0 ppm towards a Lyr (Vega, HD 172167), 8 parsecs away. Bailey et al. argue that the 
face-on debris disk around Vega is unlikely to contribute to the observed polarization signal since 
the instrument aperture is 5" and the radius of the central hole in the disk is 11" ± 2" . We show 
that the a Lyr polarization data have follow polarization trends of other nearby stars ( §5.31) so those 
data are retained in this sample. 



Santos et al.l (|201ll) surveyed interstellar polarizations towards Loop I for stars out to distances 
of 500 parsecs, using the imaging polarimeter at LNA. Typical mean accuracies lcr = 0.05% were 
achieved in the V-band. Forty stars in the Santos et al. catalog passed the regional and activity 
selection criteria above and are listed in the Appendix |B] 

Polarizations of stars used in the determination of the local ISMF direction (within 40 par- 
secs) are plotted against the distance of the star in Figure |3j For comparison, the polarizations of 
more dista nt stars (40-110 parsecs) within 90° of the heliosph ere nose are also shown (from the 
catalogs of lSantos et al.ll201ll : lHeilesll2000i : lLeroylll993Lll999|) . A general increase in polarization 
strength with distance is seen (this is discussed in more detail in Appendix 0. The polarizations 
shown in Figure [3] w ere measured i n different bandpasses. For the PlanetPol bandpass covering 
red wavelengths (|Bailey et al.ll2010|) . the Serkowski curve of polarization versus A/A max predicts 
polarizations that are up to 30% below the peak polarizations at A max (A is the wavelength of obser- 
vation and A may? whic h has a median value ~ 5500 A, is the wavelength of maximum polarization, 



Serkowski et al.lll975f) . The PlanetPol stars, located in the northern gala ctic hemisphere, includes 
the largest set of detection of very weak polarizations, while the early iTinbergenl (|1982l) dataset 
provide the most detections in the southern hemisphere. For this reason, polarization strengths are 
not used in Paper I to trace the ISMF direction. With the addition of new high- sensitivity data in 
the southern hemisphere, all polarization position angle data can be included in the study using 
weighted fits ( 34721 



Cool stars are included as target stars to provide adequate spatial coverage, but not all cool 
stars are suitable for the high- sensitivity interstellar polarizations of ~ 0.01% needed for this study. 
Magnetic activity is observed in cool stars, especially late K and M stars. With the sensitivity 
levels for most observations used in this study there is no a priori reason to omit stars with active 
chromospheres from the sample, however som e cool stars are known to show polarizations in 
chromospheric lines (IClarke & Fullertonlll996|) . Stars of spectral types G and cooler, and with 
known active chromospheres, are usually omitted from the data sample used to evaluate the local 



4 Note that the data listed for several of the Tinbergen measurements in Table 2 of iBailev et al.l d2010l) contains 
typographical errors (Philip Lucas, private communication). 
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ISMF direction. Screening for active stars was performed using the lists in lGray et al.l (|2006|) . For 
the same reason, variable stars with variations 5V > 0.06 mag are omitted to avoid other intrinsic 
magneti c phenomena, by requiring the coarse variability flag, VarFlag word H6, in the Hipparcos 
catalog (|Perryman 1997^ § to be H6=l . The analysis in Paper I did not screen for stars with known 
chromo spheric activity. 



4. Results 

4.1. ISMF direction from unweighted fit 

The first evaluation of the best-fit ISMF direction uses the same basic method as in Paper I, 
but applied to the larger set of data in this paper (©. Equation \T\ was evaluated for F{ with G n =l, 
using data points where P/dP> 2.5. All observations meeting the selection criteria are used 
in the fit, including measurements of the same star by different observers. The application of eqn. 
ED to the polarization data set (© yields an ISMF direction toward £,b=37° ± 15°, 22° ± 15° (see 
Table [T] for ecliptic coordinates). There is no significant change in the direction obtained from the 
unweighted fit when compared to the results of Paper I. The values of F{ are displayed in Figure |5] 
in both both ecliptic and galactic coordinates, and for each point in the possible ISMF grid using 
the 1° grid spacing discussed in $21 

Because of the underlying data sample, the unweighted fit emphasizes southern hemisphere 
data acquired in the 1970's. For the best-fit ISMF direction, B i=hcst , 48 stars contributed to give 
-^i=best= 0.51 for the function being minimized. Of these 48 stars, 20% have new NOT and LNA 
data and they contributed 25% of -F i=bcst , 35% have pre-1980 data from Tinbergen (1982) or Piirola 
(1977) and they contributed 29% of F i=bcst , 31% have PlanetPol data and they contributed 32% of 
-Fi=best, and 8% were observed by Santos et al. (2010) and they contributed 12% of F i= b es t- Note 
that no single set of observations appears to dominate this result. 



4.2. ISMF direction from weighted fit 

The new high-precision polarization data in the southern hemisphere, as well as additional 
high-precision data in the northern hemisphere, ($3]) makes it possible to evaluate the best-fit ISMF 
direction by utilizing all of the position angles in the sample, by weighting each value according 
by the position angle probability distribution in eqn. |2] 



5 See the VizieR catalog 1/239/. 
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If the polarizing medium were to be homogeneous over the nearest 40 parsecs in the galactic 
hemisphere (which is not the case), combining polarization data using weighted values implicitly 
weights the resulting field direction towards more distant points where the signal-to-noise (S/N) is 
larger. If in contrast there is no statistically significant field in the same spatial region, the resulting 
field direction based on weighted fits is close to the center of the spatial interval for randomly 
distributed polarization directions. 

For the weighted fits, the polarization strengths in the combined set of all data were capped at 
50 x 1CT 5 in order to minimize the possible biases due to unrecognized intrinsic stellar polarization. 
Fewer than 5% of the significant (P/dP> 2.5) data points in the fit were affected by this cap. With 
the use of weighted fits, the best-fit direction is biased towards data sets that have low quoted 
uncertainties. To avoid over-reliance on the PlanetPol data, which are quoted to a typical accuracy 
of ppm and are spatially clustered in the northern hemisphere, the P/dP value input to G n was 
capped at 3.5, and a lower limit of 1CT 5 was set on G n . With these constraints, the best-fit ISMF 
direction is oriented towards \,/3 = 260+^, 49 ± 15 or £,b=A7° ± 20°, 25° ± 20° (TableED- The 
values of Fi (eqn. [Q) are plotted in Figure |6] using G n given by eqn. [2l It is seen that the minimum 
value of Fi is more clearly defined for G Q =l than for G n given in eqn. [21 

One hundred and seventy-seven stars have contributed to this best-fit ISMF direction, 5 i= b es t- 
For this fit the mean value of eqn. [T] gives F i= b e st=485, using G n from eqn. |2] Of these 177 
stars, 26% have new NOT and LNA data and they contributed 11% of -F i= best; 39% have pre- 
1980 data from Tinbergen (1982) or Piirola (1977) and they contributed 24% of F i=best ; 12% 
have PlanetPol data and they contributed 60% of -Finest ; an d 20% were observed by Santos et al. 
(2010) and contribute 0.3% of F i= b es t- The PlanetPol contribution to the function being minimized, 
F;, is dominated by four stars that have position angles nearly orthogonal (90^^ degrees) to the 
meridians of £>i = b cs t an d that make 86% of the PlanetPol contribution. Those four stars are HD 
150680 at 10 parsecs, and three high-latitude stars at 24-31 parsecs (HD 113226, HD 116656, and 
HD 120315). 

This method of using weighted fits produces an ISMF direction that is biased by the distri- 
bution of measurement errors for diverse instruments. We have tested the effect of the weighting 
factor by raising the lower limit of G n to 10~ 4 for P/dP< 3.5. This test increases the homogeneity 
of low significance data points. Although the resulting ISMF direction is within 25° of B i=hest , it is 
closer to the center of the spatial interval indicating increased reliance on statistically insignificant 
measurements. The cap P/dP< 3.5 in G n (eqn. [2]) has the effect of deemphasizing the PlanetPol 
data where the ppm-level uncertainties otherwise dominate the result. If instead, we take G > 1CT 5 
and a Pj dP cap of 4-6, the magnetic pole is then directed towards t— 337° ± 10° and b= 28° ± 7° 
and the F minimum (eqn. [B becomes very poorly defined. We are quoting the best-fit ISMF 
direction for the weighted fit (Tabled) for values P/dP< 3.5 and G n > 1CT 5 , since these limits 
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appear to present the best overall representation of the data given the differences in measurement 
sensitivities. With future high-sensitivity data, it should be possible to treat the analysis utilizing a 
broader range of the weighting factor. 



5. Discussion 



5.1. Comparing the ISMF directions traced by the IBEX Ribbon and polarizations 

The very local ISMF direction provided by the center of the IBEX Ribbon arc is separated by 
32^3Q degrees from the direction of the best-fit ISMF (weighted fit in Table [T]). MHD heliosphere 
models that sucessfully predict the location of the IBEX Ribbon yield offsets between the center 
of the R ibbon arc and the interstellar m agnetic field that shapes the heliosphere of between 0° 



and 15° dHeerikhuisen & Pogorelov 2011 ). so the direction of the very local ISMF found here is 



marginally consistent with the direction of the ISMF that shapes the heliosphere. Polarization of 



interstellar clouds have been identified towards stars in this region of space C 


Lallement et al. 


1986; 


Frisch et al. 


2002; 


Frisch 


2003; 


Redfield & Linskv 


2008; 


Frisch et al. 


2011 


). The magnetically 



aligned dust grains causing the polarizations presumably are located in these clouds. Although 
local color excess values are too small to be measured, the variable gas-phase abundances of the 
refractory elements that make up the grains s how that several types of interstellar grains must be 



present in these clouds (|Frisch & Slavinll2012|) . The clouds are distinguished by the gas kinematics 
using the assumption of solid-body motion, and in some cases by temperature. Since most of the 
neutral ISM that is within 40 parsecs is also within 20 parsecs (@, the best-fit ISMF direction 
then may represent primarily the ISMF direction within ~ 20 parsecs of the Sun. This would 
help explain the remarkably similar directions of the ISMF traced by the IBEX Ribbon and the 
polarization data. 

If the IBEX ENA Ribbon and interstellar polarizations trace the same ISMF direction, then 
polarization strengths should increase with ENA fluxes because interstellar polarizations and the 
IBEX Ribbon both reach maximum strengths in sightlines that are perpendicular to the magnetic 
field, e.g. where B ■ R=0 for magnetic field B and radial sightline R. Figure [10] shows a plot of 
the mean value of 1.1 keV ENA fluxes that are in a 12° diameter region centered on stars that are 
within 40 parsecs, 90° of the heliosphere nose, and have P/dP> 2.5. ENA fluxes with signal- 
to-noise S/N> 3.5 are included. The data points are color-coded according to the data source in 
the left figure. Data in Figure \\0\ left, tend to be divided into four groups with different relations 
between the polarization and ENA fluxes. Group A consists of the highest quality data, such as data 
collected for this study and the PlanetPol data ($3]), and show a general increase of polarization with 
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ENA fluxes. The PlanetPol stars in group A are located in the RA = 16 H — 20 H interval. Region 
D primarily represents 1970's polarization data, where uncertainties are substantially larger. The 
trend in region B is ambiguous. The four groups of data in Figure [TOj left, are plotted in galactic 
coordinates in Figure [101 right using a color-coding that shows the data group. The increase of 
ENA fluxes with polarization for stars in group A supports the hypothesis that the ISMF sampled 
by the polarization data and Ribbon ENA fluxes are tracing the same, or similar, ISMF directions. 
Results for the remaining stars are ambiguous. 

The difference between the ISMF directions from the polarization and ENA Ribbon data may 
reveal structural components of the local ISM. The polarization data set we use to determine the 
field direction may not sample the inter stellar cloud surrou nding the Sun, the LIC. The Sun 
is within 20,000 AU of the edge of the LIC (IFrisch et al.ll201ll), and the LIC i ncludes a minimal 
amount of ISM within 90° of the heliosphere nose (e.g. Redfield & Linskyll2008 ). It is possible that 
the ISMF direction traced by the center of the IBEX Ribbon arc represents the ISMF in the LIC, 
while the ISMF traced by the polarization data represents the ISMF traced by the next cloud in the 
upwind direction, the G-cloud that is observed towards the nearest star a Cen. The LIC is ~ 22% 
ionized ( Slavin & Frischll2008l Appendix lAl). and the G-cloud is likely to be partially ionized, so 
that the ISMF will couple tightly to the gas in these clouds. The possible location of the ISMF 
traced by polarization data is discussed further below where a small rotation with distance of the 
local ISMF is found ( §5.3.21) . This rotation suggests that some of the small difference between the 
ISMF directions found from the IBEX Ribbon and polarization data may result from the volume- 
averaged ISMF direction implicit in our method for fitting the polarization data. 



5.2. Local ISMF and Loop I 



The relation between the Loop I superb ubble and t he ISM within 20 pc has been known for 
some time (see review of Frisch et al.ll201l[) . Recently, IWollebenl (|2007[) modeled the polarized 
low-frequency 1 .4 GHz and 23 GHz radio continuum of Loop I, including the North Polar Spur 
emission, as arising from two spherical shells, denoted "SI" and "S2", each with a swept up 
magnetic field traced by polarized synchrotron emission]^ An ordered magnetic field structure is 
implicit to the Wolleben model. His shell parameters suggest that the Sun is in the rim of the S 1 
shell, which is centered 78 ± 10 parsecs away towards i = 346° ± 5°, b = 3° ± 5°. Previous studies 
reached similar conclusions, although with a single superbub ble so that the center of the H° shell 
is ~ 12° south of the center of the radio continuum shell (e.g. lHeileslll998a|) . The direction of the 



6 Loop I is sometimes referred to as the 'North Polar Spur', which historically is t he brightest part of the radio 
continuum loop and juts vertically from the galactic plane near longitudes of ~ 30° (e.g jHanburv Brown et al . 1960). 



-14- 



swept-up ISMF in the S 1 shell (which has generous unc ertainties) is con sistent with the magnetic 
field direction in the local polarization patch found by iTinbergenl (|1982|) : any local fragments of 
the SI and S2 shells are not spherical, ho wever, since the Fe + column densities in local ISM do not 
trace these features defined by Wolleben ( FrischlEoiO ). The Wolleben model provides a convenient 
quantitative relation between the large-scale ISMF defined by Loop I and the ISMF close to the 
Sun. 

Local cloud kinematics, the local ISMF, and local gas-phase abundances all suggest an origin 
related to Loop I. In the local standard of rest (LSR), the local cluster of clouds flows p ast the Sun at 
a bulk velocity of V — —17 km s _1 , and towards us from the direction 335°, —5° (|Frisch et al. 
20111). Solid-body motion is assumed in determining this vector. The flow direction is ~ 14° from 
the center of the SI shell, which is located at 335°, —5°. Although the SI shell is an evolved 
superbubble remnant, the local fragment of the shell appears to still be radially expanding into the 
low density region of the Local Bubble. The best ISMF direction from the weighted fit (Table 1) 
makes an angle of 76° with the local ISM bulk flow vector. These numbers are consistent with a 
scenario where the normal to the shell rim is parallel to the ISM velocity, and the magnetic field 
is compressed in the shell rim so that the field direction is nearly perpendicular to the normal of 
the shell surface. An association of local ISM with a supernova shell is consistent with the high 
and variable abundances of the refractory elements Fe, Mg, and Si in local ISM, which sugges t 
inhomogeneous grain destruction by interstellar shocks (|Frisch et al.ll2011uFrisch & Slavinll2012|) . 



5.3. Ordered ISMF near Sun 



A detailed study of the behavior of interstellar polarizations and ISMF behavior with distance 
is possible with the sensitive PlanetPol dat a for stars located in the interval RA = 16 H — 20 H , 
which contains the RA> 17 H region where bailey et~aL (|201o|) found an increase of polarization 
with distance. This region overlaps the low latitude part of the North Polar Spur, where aligned 
optical polarization vect ors define the ISMF for distances 60-120 parsecs (|Mathewson & Ford 
1970j; ISantos et al.ll2011r) . The behavior of the ISMF for distances < 60 parsecs in the North Polar 
Spur direction is provided by PlanetPol data. 



5.3.1. Polarization versus distance and column density in 16 H — 20 H region 



Bailey et al.1 (120 10f) found that the increase of polarization with distance was greater for stars 



with RA> 17 H than for stars in other spatial intervals. The increase of polarization with star 
distance for PlanetPol stars in the RA = 16 H — 20 H region and with P/dP> 3.0 is shown in 
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Figure|7] right. For this set of uniform quality data a maximum value for the polarization is found 
for each distance interval within ~ 60 pc, and there are no observed points that lie above the 
upper envelope marked by the filled points. Polarizations of stars beyond that follow this trend. 
Stars which defin e this maximum are often refereed to as the polarization 'upper envelope' (e.g. 



Anders sonll2012l also see Appendix 0, and we use that term here. In general, polarization data 
tend to have an upper envelope when plotted against either distance or color excess due to either 
depolarization or patchy ISM. A rough estimate of the dependence of polarization on distance 
can be found by first binning stars into 10 parsec intervals, and then selecting the stars with the 
maximum polarization in each interval. Omitting bins with no stars, seven stars (filled circles in 
Figure \7} can be used to define the upper limit of polarization as a function of distance. A linear fit 
though the envelope stars gives the approximate distance dependence of the maximum polarization 
in the RA = 16 H — 20 H region of 

logP ~ 1.200 + 0.013 * D (3) 

where P is polarization in ppm and D is the distance in parsecs. The closest envelope star is a 
Lyr at 8 parsecs (Table [3] lists the envelope stars). The locations of the RA = 16 H — 20 H stars 
are plotted in galactic coordinates in Figure [7] left. Two unpolarized PlanetPol stars in this region 
are not plotted, HD 163588 and HD 188119, both of which are located at larger longitudes than 
the other stars. Stars that form the upper envelope of the polarization-distance relation are loosely 
confined between galactic longitudes of 30° and 80°, and latitudes of -5° to 30°. 

The upper envelope of the polarization-distance relation (Figure[7|) can be used to estimate the 
gas density gradient along the ISM holding the ISMF, since dust and gas have similar spatial distri- 
butions in diffuse clouds. Few of the PlanetPol stars have measured H° or D° column densities (e.g. 



Wood et al.ll2005|) . so column densities are estimated by converting polarization to color excess, 
and then converting color excess to A r (H), using standard relations for diffuse clouds. The polar- 
ization for dust grains aligned by an ordered ISMF, with no rando m component, provi des an upper 



envelope of the P vs. E(B-V) relation of P(%) = 9E(B - V) dFosalba et al.ll2002h . Interstellar 



gas and dust have similar spatial distributions in the low opacity ISM. Copernicus observations 
of the saturated H° Lya absorption line and molecu lar hydrogen lines showed that N(H)/E(B- 
V)=A^(H°+ H 2 )/£(B-V)= 5.8 x 10 21 - cm" 2 mag" 1 (bohlin et al.lll978h FI The Copernicus value 



underestimates hydrogen column densities for partially ionized gas since independent measure- 
ments of ionized gas are not available for most sightlines. The relation between polarization and 
column density is estimated by combining the distance-polarization fit with the empirical P vs. 



7 We use the ratio for the ISM towards all stars, rather than the ratio for the interc loud medium, 7V(H°+ H2)= 
5.0 x 10 21 -E(B-V) cm" 2 , that might be relevant to the nearest ISM feohlin et al.lll978b . 
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E(B-Y) and E(B-V) vs. 7V(H) relations, giving 

logiV(H) = 18.01 + 0.013 * D. 



(4) 



for stars tracing the upper envelope of the P-E(B-V) relation and P-distance relation. For other 
stars, foreground depolarization can lead to lower values for iV(H ). This gas-distance relation can 
be c hecked aga i nst th ree nearby stars in the RA = 16 H — 20 H region with iV(H°) values determined 
bv lWood et all t005\) . HP 155886 (36 Oph), 165341 (70 Oph), and HD 165185. The ratio of the 
predicted iV(H) (eqn. |U) to the observed iV(H°) for these three stars is 1.35 ± 0.35, indicating 
depolarization and/or inhomogeneous ISM is present nearby. 



5.3.2. Polarization versus position angle in 16 H — 20 H region 
Optical polarization data that traces the H° filamentary structure of Loop I beyond ~ 8 



parsecs reveal an ordered ISMF that follows the axis of an H I filament of Loop I (IHeiles 1998a). 
Optical polarization vectors (Figure HJ, and the location of the RA = 16 H — 20 H region towards 
the North Polar Spur, suggest that the very local ISMF direction may also be ordered. Figure 
[8] displays the polarization position angles, PA R a, for PlanetPol stars with P/dP> 3.0 in this 
interval, as a function of distance. With the exception of the outlier star HD 164058, which has 
PAra=145° ± 0.5° and a distance of 45 parsecs, the envelope stars within 55 parsecs display 
a smooth rotation of polarization position angle with distance. HD 164058 is a K5III star and 
belongs to a spectral class where intrinsic polarization occurs because of chromospheric magnetic 
activity and was omitted from the fit in eqn. [5] as an obvious outlier (see below). A linear fit to 
position angles versus distance gives 

PA RA (D) = 35.97(±1.4) - 0.25(±0.03) • D (5) 

and a x 2 of 0.559, where PA R a is m degrees and the distance D is in parsecs. The small value of 
X 2 per degree of freedom (2) suggests that the uncertainties on PAra quoted by lBailey et alj(l2010|) 



may be generous. @ The probability that the computed x 2 value would have a value larger than 
0.559 is 0.76, indicating that parameters of the fit in eqn. [5] are believable. After excluding HD 
164058, the relation between the polarization direction and star distances of the four envelope stars 
within 55 parsecs appears to be linear. The large difference of 7.8cr pp between the position angle of 
the outlier star and the best-fit line (eqn. |5|), where a pp is the mean difference between the best-fit 



iBailev et alj (120101) include several different effects when calculating uncertainties for their measurements, includ- 
ing counting statistics, possible intrinsic telescope polarization, and in some cases atmospheric polarization caused by 
Sahara dust. 
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line and the position angles of PlanetPol stars within 55 pc except for the outlier, justifies omitting 
HD 164058 from the calcuation of eqn. [5j The envelope stars are partly self- selected by the fact 
they represent sightlines where intervening depolarization is minimized, so the linear fit to the 
polarization position angles suggests the existence of a nearby ordered ISMF that rotates ~ 0.25° 
perparsec. This r otation is confined to the small region of space (filled dots in Figure 7, left) where 
Baile y et al. found the systematic increase in polarization strengths with distance (Figure [7] 

right). The envelope stars appear to sample the uniform part of an extended medium that contains 
a magnetic field and reaches to within 8 pc of the Sun. 

The standard deviation of position angles around the ordered component is a measure of the 
magnetic turbulence in this region. For all stars in the RA = 16 H — 20 H region and within 55 
parsecs, excluding the outlier star, the standard deviation of polarization position angles about the 
best-fit line is 23°, which we interpret as magnetic turbulence. 

The smooth variation of the ISMF direction over distance suggests that the four nearest enve- 
lope stars are embedded in the polarizing medium, which extends between a Lyr (Vega, 8 parsecs) 
and 5 Cyg (52 parsecs). The sightline s to all of the env elope stars sample either the SI shell, the S2 
shell , or both sh ells, according to the Iwollebenl (|2007|) model and the shell configurations plotted 
in lFrischl (|2010|) . The angular spread on the sky of the envelope stars is 50°, which suggests that the 



ISMF in the nominal SI and S2 shells consists of a homogeneous "matrix" component of smooth 
material that is sampled by the envelope stars, and patchy material and/or a turbulent ISMF that 
produces depolarization of, and turbulence in, the position angles of non-envelope stars. 

If the polarizations of the envelope stars are tracing matrix ISM in this region, this mate- 
rial should be detectable as interstellar absorption features. Absorption lines have been measured 
towards the four envelope stars in th e fit. All fou r stars show absorption at or close to the ve- 
locity of the Apex Cloud defined by Frisch ( 2003|) F1 so it is tempting to assume that the 'Apex 
cloud' represents the smooth ISM filling the local region of the SI and S2 shells. Interstellar 
absorption lines have also been measured towards a fifth polarized star, HD 187642, in region 
RA = 16 H — 20 H , where a component at the Apex cloud velocity is also found. Towards the 
envelope star HD 159561, the Apex cloud component is characterized by weak Ca + absorption 
and a typical ratio Ti + /Ca + =0 .41, suggesting it contain s some neutral gas since Ti + and H° have 
the same ionization potential (IWeltv & Crowtherll2010l). A mor e recent and complete cloud iden- 
tification scheme was developed by iRedfield & Linskvl (|2008|) . where three of the Apex cloud 
components towards envelope stars are instead assigned to the MIC cloud , while the fourth com- 
ponent is left as unassigned. Three of the four stars in IRedfield & Linskvl (I2004bh that are in the 
RA = 16 H — 20 H region have components at the Apex velocity, however the cloud temperature 



9 This cloud vector is similar to the "Panoramix" cloud that was identified much earlier by 



Lallement et al 



(1986). 
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ranges from ~ 1, 700 — 12, 500 K for these components. A second possible carrier of the ordered 
field is the 'G' cloud since it has a wide angular extent in the upwind direction, is within 1 .3 parsecs 
of the Sun, and is observed towards the envelope stars. The identification of the cloud carrying the 
ordered component of the ISMF requires additional study. 



5.3.3. Polarization versus interstellar radiation field in 16 H — 20 H region 



Radiative torques are known to play a role in th e alignment of dust grains in dense clouds 
where the radiation field has a preferred direction (e.g. lAndersson & Potterll200~6 ; Draine & Weingartnen 



19971 : lAndersson et al.ll201 it Lazarian & Hoangll2008|) . but this does not appear to be the case lo- 
cally. The interstellar radiation field extrapolated to the location of each star in the RA = 16 H — 20 H 
region appears to anticorrelate with polarization strengths. We believe that this correlation is a 
chance coincidence of the relative locations of the brightest radiation sources and polarized stars. 
Figure [9] shows polarizations for stars inside (filled symbols) and outside (open symbols) of the 
RA = 16 H — 20 H region, plotted against the total far-UV 975A flux that has been extrapolated 
to the position of each polarized star, while ignoring any possib le opacity effects. T he 975A ra- 
diation flux is given by the 25 brightest stars at this wavelength (|Opal & Wellerill984|) . which are 
primarily located in the third and fourth galactic quadrants because of the geometry and low opac- 
ity of the Local Bubble ( see, e.g., the interstellar radiation field at the Sun from the S2 UV survey, 
Gondhalekar et al.lll980|) . These polarization data show a trend for linear polarization to increase 
as radiation flux decreases. The rough anticorrelation between radiation fluxes and polarizations 
is shown by the fit log P= 4.61 — 0.054*F 3 (dashed line in Figure [9]), where P is in units of 
10~ 6 and F 3 is the radiation flux at 975 A in units of 10 3 photons cm~ 2 s _1 A -1 . This apparent 
anticorrelation appears to be the result of the relative locations of the brightest hot UV-bright stars 
and the most distant polarized stars, that define the low-flux region of the trend (see Appendix |TJ 
for more details). 



5.3.4. CMB dipole anisotropy and the local ISMF 



Large scale global deviations from a uniform blackbody spectrum of the CMB ra diation define 
a dipo le anisotropy that is oriented towards £= 263.85° ±0.1° and b= 48.25° ±0.04° dBennett et al. 



20031) . The great circle on the sky that is midway between the two poles of the CMB dipole 
displays a striking symmetry around the directi on of the hel iosphere nose direction defined by the 
flow of interstellar He° through the heliosphere ( Frischll2007 ). The best direction of the heliosphere 
nose is defined by the new IBEX consensus directi on for the interstellar H e° flow, which includes 
correction for propagation and ionization effects (|McComas et al.ll2012l . see Appendix lAl). This 
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new consensus IBEX direction for the heliosphere nose more closely aligns the nose direction with 
the great circle that is midway between the hot and cold poles of the CMB dipole moment. The 
angle between the heliosphere nose and dipole anisotropy direction is 91.5° ±0.6°. The probability 
that the great circle that is evenly spaced between the poles of the CMB dipole anisotropy, which 
is defined to within ~ ±0.10°, would pass through the ~ 3° x 3° area on the sky that includes the 
heliosphere nose direction, is less than 1/200. Figure [TT| shows the great circle in the sky that is 
equidistant from the antipodes of the CMB dipole moment (gray line), as well as the hot (red) and 
cold (blue) dipole directions. 

The 90-degree great circle midway between the dipole directions also passes through the 
best-fit ISMF direction, when uncertainties are included (Table [JJ). The ISMF direction from the 
weighted fit us plotted as a purple dot in Figure [TT] together with the uncertainties on this direction 
(purple polygon) and the 90-degree circle. That the great-circle dividing the antipodes of the 
CMB dipole moment would also pass through the ISMF pole, which is defined by an area of 
~ 300 square-degrees, including uncertainties, is less than one in 10 4 . A deeper knowledge of 
the configuration and properties of the global heliosphere and local ISMF will elucidate whether 
additional CMB foregrounds might be present. 

Although a possible relation between the CMB dipole and heliosphere may seem far-fetched, 
a number of studies have shown that the symmetries traced by the low-£ CMB moments are not 
consistent with the predictions of gaussianity of the universe based on standard cosmological mod- 
els, and in particular the quadrupole and octopol e moments show evidence of the ecliptic geometry 
(e.g. ISchwarz et al.ll2004t IStarkman et al.ll2009l and citations in and to these articles). Since the 
ISM flows through the heliosphere from a direction 5° above the ecliptic plane, these low-£ sym- 
metries are equivalent to symmetries around the heliosphere nose (e.g. see Figures 6 and 7 in 
Frischll2007l) . Unrecognized CMB foregrounds related to the heliosphere may be present, with the 
most likely candidate being emission from i nterstellar dust inside the heliosphere since the dust 
traces both the gas flow and ISMF direction dsiavin et al.ll2012l : lFrisch & Slavin|[2012l:lsiavin et al 
20091) . 



5.4. Speculative implications 



5.4.1. Local ISMF: Interarm or Arm ? 



Comparisons between the ISMF at the heliosphere and the ISMF observed in other nearby 
regions provides insight into the uniform and random components of the local ISMF, such as 
whether the field has the large-scale coherence of interarm fields or of the small-scale random 
fields seen in spiral arm regions. Random components of the ISMF in spiral arms are coherent over 
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scale sizes of 5-50 parsecs; Faraday rotation structure funct i ons indicates a more tu rbulent ISMF 
with outer scales of 10-20 parsecs (|Rand & Kulkarnilll989t lHaverkorn et al.ll2006|) . In contrast, 
Faraday rotation measures for interarm sightlines show a more uniform component that is coherent 
over scales of 100-200 parsecs. The ratio of the uniform to random components of interstellar 
polarizations in the p lane of the sky is typically ~ 0.8, when global polarization data are considered 
iFosalba et ap2002h . 



If the local ISMF sampled by the IBEX Ribbon and optical polarization data is part of an ISMF 
that is coherent over several hundred parsecs, the coherence should appear in the ISMF direction 
obtained from Faraday rotation measure (RM) data on nearby p ulsars that yiel d an electron-density 
weighted measurement of the parallel component of the ISMF. Salvati ( 201o|) fit the RM and dis- 
persion measures of four nearby (distances 160-290 parsecs) pulsars that sample the low density 
interior of the Local Bubble in the third galactic quadrant (£~ 180° - 270°). The ISMF strength 
obtained from these data is 3.3 /xG, and the field is directed towards 5°, 42°. The third galactic 
quadrant is opposite to the RA = 16 H — 20 H region in the sky. The difference in the direction of 
the ISMF traced by the RM data and IBEX Ribbon, ~ 23°, is easily accommodated by the field 
curvature indicated by the envelope stars in the RA = 16 H — 20 H region (including a Lyr 8 parsecs 
away), and may indicate that the ISMF around the Sun extends into the third galactic quadrant. 

The strength of the ISMF that shapes the heliosphere, ~ 3 fiG, is close to the field strength 
found from the pulsar data. At the heliosphere, magnetic field strength is estimated from the 
pressure equilibrium between the inner heliosheath ions traced by ENAs and the ISM, and from 
the magnetic tension require d to produce the 44° off set between the observed ENA minimum and 
downwind gas flow direction ( Schwadron et al.lboill). while in the LIC it is found from equilibrium 
between gas and magnetic pressures (|Slavin & Frischll2008|) . The polarity of the ISMF found from 
the pulsar RM data is directed towards A,/3=232°,18°, and therefore is directed upwards through 
the ecliptic plane (see Table 1 for galactic coordinates). Although the field direction is consistent 
with the local ISMF shaping the heliosphere, the polarity appears to be the opposite of the polarity 
built into most MHD heliosphere models (Appendix |A"1). 

In the opposite direction of the sky from the pulsars, the RA = 16 H — 20 H region tests the 
coherence of the ISMF along the sightlines that are tangent to the Loop I shell. The slow rotation of 
the nearby ordered field component and its extension to within 8 parsecs pf the Sun, together with 
the similarity of ISMF directions from the IBEX Ribbon and pulsar data, suggests that the Loop I 
ISMF may extend into the third galactic quadrant. Figure [T]shows that the dust cavity caused by the 
Loop I expansion covers the northern polar cap and extends into the third galactic quadrant at high 
latitudes. The coherence of the ISMF can be found by comparing values from these three regions: 
(1) at the solar location, from the center of the IBEX Ribbon arc (to within ~ 15°, see Appendix 
|A"1); (2) the nearest 40 parsecs in the galactic center region, including the tangential direction of 
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Loop I where the ordered component of the ISMF is found; (3) the low density void of the Local 
Bubble in the third galactic quadrant. The Sun is located (by definition) in region 1, and between 
regions 2 and 3. The difference in the ISMF directions determined from regions 1 and 2 is ~ 32°, 
and it is ~ 22° between regions 1 and 3. The field curvature must extend from region 2, through 
the solar location in region 1, and into the third galactic quadrant corresponding to region 3. These 
data suggest that the local ISMF at the heliosphere is part of a large scale field that extends several 
hundred parsecs in space, such as is characteristics of interarm regions. The Loop I ISMF may be 
"opening out" into an interarm field as it expands into the low density third galactic quadrant. 



5.4.2. Galactic cosmic rays and the local ISMF 



Galactic cosmic rays (CR) in the energy range of 50 GeV - 40 TeV trace the local ISMF 
because the gyroradius is on the order of several hundred AU in a 1 fiG magnetic field. Large 
scale anisotropics an d smaller scale excess regions, have been detected in the distribution of CRs 



in this energy range ( Nagashima et aL 



998; Hal l et al.lll999blAbdo et al. 2008: Abbasi etall20U 



Surdo & Argo-Ybj Collaborationll2011|) . It has been suggested that the Sun is located in a magnetic 
flux tube in th e ISM, and that such a flu x tube would be a conduit for galactic cosmic rays (GCRs) 
(IFrischl Il997|) . ICox & Heleniusl (|2003|) presented a model of the Local Bubble that placed the 
Sun in a flux tube that has detached from the walls of the Local Bubble. iSchwadron et al.l (|2012|) 
have identified anisotropics in the CR distribution towards the heliotail direction caused by the 
propagation of CR's through the magnetic structure surrounding the heliosphere. Therefore it is an 
interesting exercise to compare the GCR asymmetries with the local ISMF direction. 



Using cosmic ray fluxes in the northern and southern hemispheres, Nagashima et al.l (11998) 
described two anisotropics in terms of abroad sidereal galactic anisotropy, and a narrow anisotropy, 
the "tail-in" anisotropy, described by a maximum towards the anti-apex direction of the solar mo- 
tion0 The galactic anisotropy arrives from the direction a, 5 = H , —20° (£, b ~ 61°, —76°). The 
tail-in excess corresponds to the direction a, 8 — 6 H , —24° (£, b ~ 230°, —21°), with a half-width 
of ~ 68° and a maximum effect near ~ 10 3 GeV. The tail-in excess was originally attributed to 
the h eliotail because i t is approximately opposite to the direction of the solar motion through the 
LSR. lHall et alJ (119991) fit the anisotropics as Gaussian features, and found a direction of maximum 
anisotropy for the tail-in excess towards a, 5 ~ 4.6 H , —14° (I, b ~ 210°, —36°), and a loss-cone 
centered at a, S ~ 13 H , 0° (£, b ~ 123.°, 63.°). To within the uncertainties, the direction of the 



10 The solar apex motion corresponds to a velocity of 18 ± 0.9 km s _1 directed towards I ~ 47.9° ± 3.0°, b = 
23.8° ± 2.0°, however these values changes on a fairly regular basis with improved understanding of the astrometric 
data for the comparison stars used to establish the solar motion. 
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tail-in excess found by Hall et al. overlaps both the downwind He° flow vector that defines the tail 
of neutral ISM gas behind the heliosphere, and the axis of the ISMF direction determined from the 
weighted fit (towards i,b= 229° ± 30°, -28° ± 30°). Figure H] shows the directions of the ISMF, 
tail-in excess compared to the ISMF direction, and the true h eliotail corresponding to the ENA flux 
minimum centered ~ 44° west of the downwind gas flow (|Schwadron et al.ll201 1|) . 



6. Summary and conclusions 

We have determined the direction of the interstellar magnetic field within 40 parsecs of the 
Sun and 90° of the heliosphere nose (~ 15° from the galactic center) using starlight polarized 
by magnetically aligned interstellar dust grains. The analysis is based on the assumption that the 
polarization E-vecto r is parallel to the direction of the ISMF and builds on the analysis in Paper I 



(|Frisch et al.H2010al) . The new technique for finding the ISMF direction minimizes the mean of an 
ensemble of sines of polarization position angles ($2]). The use of an expanded data set, including 
new polarization data that have been acquired at several observatories for this study (@, permits 
the use of weighted fits to the ISMF ($4]) as opposed to the unweighted fit used in Paper I. The 
best-fit ISMF from the weighted fit is towards £,b= 47° ± 20°, 25° ± 20°. The direction obtained 
from the unweighted fit is consistent with the results of the weighted fit (Table [T]). 

The best-fit ISMF directions obtained from the weighted fit and from the center of the IBEX 
Ribbon arc (Table 1) are separated by 32+^ degrees ( ^5.1| ). Although heliosphere models indicate 
that there could be a small offset between the direction of the ISMF shaping the heliosphere and 
the Ribbon arc center (0 — 15°), the agreement between these two field directions is remarkable 
considering that the polarization-based measurement is determined from light of stars up to 40 
pc away. The similarity of these directions encourages the conclusion that the nearby ISMF is 
coherent over decades of parsecs, while the difference indicates that the local ISMF either has a 
curvature or is turbulent. 

Interstellar optical polarizations reach maximum strengths in directions perpendicular to the 
ISMF, and the IBEX Ribbon appears in directions perpendicular to the ISMF draping over the 
heliosphere. Hence if the directions of the best-fit ISMF and magnetic field traced by the IBEX 
Ribbon are related, the polarization strengths and ENA fluxes will be proportional. In one region 
of the sky that touches the Ribbon, where the ISMF is observed within 8 pc of the Sun, ENA fluxes 
and polarization strengths are found to increase together ( ^5.11) . This effect is seen only in the most 
sensitive polarization data. 



In one region of the sky in the right-ascension interval RA = 16 H — 20 H . lBailey et al.1 (I2010|) 



used their high-sensitivity PlanetPol data to conclude that polarization systematically increases 
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with distance ( §5.3.11) . Using stars that form the upper envelope of the polarization versus distance 
relation in that region, we find that logP = 1.200 + 0.013 * D. For standard relations between 
polarization and color excess, and color excess versus iV(H), the distant-dependent increase in col- 
umn density in this spatial interval becomes log A(H)= 1 8. 01 +0.01 3 *D; the relation overpredicts 
the actual iV(H°) measured for several stars by ~ 35%. 

This same region, which overlaps a low-latitude portion of Loop I, appears to contain an 
ordered ISMF that extends to within 8 parsecs of the Sun ( §5.3.2| ). Stars forming the upper enve- 
lope of the polarization-distance relation define a unique group where foreground depolarization 
is minimized. For the nearby envelope stars, the polarization angle depends on the star distance 
with the relation PA RA = 35.97(±1.4) - 0.25(±0.03) ■ D The ISMF rotates by ~ 0.25° 

pc _1 , yielding a 10° ISMF rotation over a 40 parsecs thick magnetic layer. The nearest polarized 
star showing the ordered ISMF is HD 172167 (a Lyr, 8 parsecs), where the polarization direction 
is consistent with the direction of best-fit ISMF field. The dispersion of position angles about this 
nearby ordered ISMF component suggests that the turbulent ISMF component is ~ ±23°. These 
results are conditioned on the exclusion of a 7.8cr outlier star from the analysis, and on the small 
number of stars (4) available for defining the position angle rotation. Confirmation of this result 
would require observations at ppm accuracy towards very faint stars, V> 6 mag, in the spatial 
interval £~ 25° - 70°, 6~ 0° - 325°. 

The best-fit magnetic field direction, the flow of local ISM past the Sun, and variable abun- 
dances of refractory elements in the local ISM all suggest that the local ISM is associated with an 
expanding fragment of the SI shell of the Loop I superbubble ( §5 .21) . The bulk LSR flow of the 
local ISM past the Sun has an upwind direction within 15° of the center of the S 1 shell. The best-fit 
ISMF makes an angle of ps 76° with the flow direction, suggesting the field is parallel to the rim 
and perpendicular to the flow velocity. Such a configuration would be expected for an ISMF swept 
up and compressed in an expanding superbubble shell. Variable abundances of refractory elements 
in local ISM indicate recent shock destruction of some of the local dust grains. A consistent inter- 
pretation of the kinematics of the local ISM and the direction of the best-fit ISMF would be that the 
ordered component represents the ISMF that was been swept up in the expanding SI superbubble 
shell. The flow of ISM past the Sun is deceleratin g, as shown by the deviations of velocity com- 



ponents from a rigid-body flow (|Frisch et al.ll201 1|) . Cloud collisions in the decelerating flow may 



generate the observed magnetic turbulence of ±23°, with the ordered component representing the 
matrix ISM between the clouds. 

Polarization strengths for the envelope stars in region RA = 16 H — 20 H anticorrelate with 
the interstellar radiation field at 975 A ( §5.3.31 ). This anticorrelation is probably due to the fact that 
those stars are in the first galactic quadrant and, by chance, the most distant of those stars are also 
the most distant from the bright stellar far ultraviolet and extreme ultraviolet radiation sources that 
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are mainly located in the third galactic quadrant. 

The ISMF direction from pulsars in the third galactic quadrant is within 22° of the ISMF 
direction from the IBEX Ribbon. A common field strength of ~ 3 /iG is found from the pulsar data, 
a range of heliosphere diagnostics including plasma pressure and the deflection of the heliotail, and 
models of the ionization of the ISM surrounding the heliosphere. This suggests that the ISMF near 
the Sun is coherent over large spatial scales, such as expected for interarm regions. One source 
of uncertainty is the polarity of the field, since the polarity derived from pulsar data is opposite in 
direction from the polarity that is assumed by several of the MHD heliosphere models. 

The new ISMF direction, together with the new consensus direction for the ISM flow through 
the heliosphere derived from IBEX observations of interstellar He°, strengthens the spatial coinci- 
dence between the geometry of the CMB dipole moment, the heliosphere nose, and the local ISMF 
C §5.3.4b . The great circle that divides the two poles of the CMB dipole moment passes within 
1.5° ± 0.6° of the heliosphere nose defined by the IBEX measurements of the flow of interstel- 
lar He° in the heliosphere, and passes through the best-fit ISMF determined by weighted fits to 
polarization position angles, to within the uncertainties. 

The axis of the best-fit ISMF direction from polarization data extends through the direction of 
the 'tail-in' excess of GeV-TeV cosmic rays that has been attributed to the heliotail. The direction 
of the tail-in excess coincides poorly with the flow of interstellar He° through the heliosphere. 

High- sensitivity interstellar polarization data provide an opportunity to explore the ISMF and 
magnetic turbulence in the very diffuse ISM close to the Sun, and our very local galactic environ- 
ment. The implications of understanding this field may affect our understanding of the most distant 
reaches of the Universe. The domain of the heliosphere and the local ISM are filters through which 
all observations about our more distant universe must pass. 
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A. Appendix: Heliosphere models and the IBEX Ribbon 



Helio sphere models predict the density and ionization of the LIC and the ISMF that shapes the 



helio sphere (e.g. lOpher et al.l 120091: IPogorelov et al.ll2009at iRatkiewicz et al.l 120081 : IPrested et al 



20101), and originally showed that the Ribbon forms where the sightline is perpendicular to the di- 
rection of the interstellar field draping over the heliosphere ( Schwadron et al]|2009[ IHeerikhuisen etal 
2010 ; Ichalov et allboid : iGrvgorczuk et al.lboi ll iHeerikhuisen & Pogorelovlboi ll : IRatkiewicz et al 
2012|) . The exact mechanism generating the Ribbon is a c onundrum, partly because the properties 



of tu rbulence upstream of the heliopause are unknown (|Gamayunov et al.1 1201 Oc iFlorinski et al. 



201 0), so that several possible scenarios for the formation of the Ribb on have been suggested (e.g. 
McComas et al.ll2010l : ISchwadron et al.ll201ll : lGrzedzielski et al.ll2010|) . One possible Ribbon for- 
mation mechanism requires outflowing ENAs to escape from the heliosphere, and charge exchange 
with interstellar protons in the outer heliosheath. These ENAs subsequently create secondary in- 
flowing ENAs from another quick charge exchange with interstellar neutrals before the ring-beam 
distribution of the ions (with a pitch angle of ~ 90°) i s disrupted by turbulence, and t herefore 
giving the B ■ R— alignment. Testing for such a model. IHeerikhuisen & Pogorelov ( 201 lb varied 
the strength and direction of the ISMF, and the interstellar neutral H and proton densities, to ob- 
tain a best-fit to the geometry of the Ribbon. These models suggest the Ribbon is formed roughly 
~ 100 AU beyond the heliopause, by a local ISMF that is directed away from the galactic coor- 
dinates £,6=33° ± 4°, 53° ± 2°, for a field strength of 2-3 fiG. For comparison, an earlier MHD 



This preprint was prepared with the AAS L^TgX macros v5.2. 
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model of the heliosphere asymmetries implied by the deflection of the solar wind in the inner he- 
lios heath gave an ISM F directio n of £= 10° - 22°, b= 28° - 38°, and field strength ~ 3.7 - 5.5 



fiG (lOpher et al.ll2009|) . For the iHeerikhuisen & Pogorelovl (|201 11) models, the angle between the 
center of the Ribbon arc and the ISMF direction is ~ 0° — 15°. 

The global heliosphere provides an in situ diagnostic of the magnetic field and plasma shaping 
the heliosphere. The ~ 48° offset between ISMF direction given by the Ribbon arc center and the 
velocity of the inflo wing neutral interstellar gas creates asymmetries that can be predicted using 



MHD models (e.g. lPogorelov et a l. 2009bl: lQpher et alJ l2009: Rat kiewicz et al.ll2008l : IPrested et al. 



2010|) . The data tracing the heliosphere asymmetries include the different distances of the solar 



wind termination shock found by Voyager 1 i n the northern ecl i ptic hemisphere and that f ound 



by Voyager 2 in the south (94 AU vs. 84 AU, IStone et all 12008 



Richardson & Sto ne 2009&, the 



~ 8° offset between the upwi nd directions of the neutral interstellar H° and He° flowing into the 
heliosphere (using data from Lallement et al.l 12010c. iMcComas et al.ll2012l) . and the preferential 
alignment exhibited by the 3 kH z plasma emissions from the ISM beyond the heliopause that were 
detected by both Voyagers (e.g. iKurth & Gurnettll2003|) . The ISMF directions from these models 
are generally consistent with the field direction from the IBEX ribbon. 



McComas et al.l (120 12|) have determined a new consensus direction for the velocity vector of 
interstellar He° flowing through the heliosphere. The He° gas velocity, 23.2 ± 0.3 km s -1 (Table 
CO, agrees with the (less precise) velocity of interstellar dust traveling through the heliosphere, 
after correcti on for propagation effects, of 24.51H km s _1 from galactic coordinates i,b= 8° ± 
15°, 14° ±4° dFrisch et all 1 9991 : iKimura et al.ll2003l) . It also agrees with the velocity of interstellar 
H° interacting with the heliosphere found from the first spectrum of Lya emission from interstellar 
H° inside of the heliosphere, 22.5 ±2.8 km s" 1 , after c orrection for H° propagation and to the 
consensus IBEX upwind direction (|Adams & Frischll 1 977[) . 



The outer boundary conditions of the heliosphere models are set by the physical properties 
of the LIC, which must be reconstructed from ISM observations using photoionization models be- 
cause of the low opacity of the cloud. Ribbon models are highly sensitive to the interstellar bound- 
ary conditions of the heliosphere; 15% variations in the interstellar param eters lead to pronounced 
differences in the location and width of the Ribbon (|Frisch et al.ll2010b|) . Photoionization models 
indicate that the LIC is a low de nsity, n ~ 0.27 cm" 3 , partially ionized cloud, H + /H ~ 22% and 
He+/He ~ 39% (Model 26 in Isiavin & Frischll2Q08b . Observations of ISM inside of the helio- 



sphere also probe the LIC properties. When the IBEX measurements of interstellar 0°, and Ne° are 
corrected for propagation effects through and entering the heliosphere, and ionization in the LIC 
using photoionization models, the interstellar abundance ratio 0/Ne= 0.27±0.10 that is recovered 
agrees with other local ISM measurements and suggests that over 50% of the O is depleted onto 
interstellar dust grains (|Bochsler et al.ll2012|) . IBEX and Ulysses measurements of interstellar He° 
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indicate the LIC is warm, 6300 ± 390 K dWittdl2004l : McComas et alJl2012h 



New IBEX-LO observations of the flow of interstellar He° through the inner heliosphere 



(IMcComas et al.1 2012|) indicate that ISM flow is ~ 12% slower than found from earlier Ulysses 
data dwittel 2004), which lowers the interstellar ram pressure on the heliosphere by ~ 22% and 



may affect the models of the Ribbon location. The value of the He° flow vector is given in Table 
1, together with the la uncertainties for longitude, latitude, and velocity. The new IBEX upwind 
direction has shifted to 3.6° east of the Ulysses direction. The uncertainties are somewhat larger 
when the bounding range constrains are included; these constraints set the limits as £= 5.25l^ 
degrees, b= 12.03jjjg degrees, the velocity to — 23.20^J km s _1 , and the temperature to 5000- 
8300 K. 



B. Appendix: List of stars from the literature included in study 

The stars within 40 parsecs from the PlanetPol catalog included in this study are: HD 97603 
(HR 4357), HD 140573 (HR 5854), HD 139006 (HR 5793), HD 116656 (HR 5054), HD 156164 
(HR 6410), HD 95418 (HR 4295), HD 112185 (HR 4905), HD 161096 (HR 6603), HD 177724 
(HR 7235), HD 127762 (HR 5435), HD 153210 (HR 6299), HD 120315 (HR 5191), HD 113226 
(HR 4932), HD 1 12413 (HR 4915), HD 163588 (HR 6688), HD 95689 (HR 4301) and HD 131873 
(HR 5563). 



The nearby stars from iSantos et al.l (|2011[) catalog that have been used in this study include: 
HD 223889 (HIP 117828), HD 120467 (HIP 67487), HD 173818 (HIP 92200), HIP 65520, HD 
127339 (HIP 70956), HIP 82283, HIP 50808, HIP 106803, HD 162283 (HIP 87322), HIP 1 14859 
HIP 83405 HIP 96710 HIP 90035 HIP 87745 HD 176986 (HIP 93540), HD 175726 (HIP 92984), 
HIP 61872, HIP 95417, HD 155802 (HIP 84303), HD 166184 (HIP 88961), HD 161098 (HIP 
86765), HD 182085 (HIP 95299), HD 164651 (HIP 88324), HD 119638 (HIP 67069), HD 125184 
(HIP 69881), HD 95521 (HIP 53837), HIP 66918, HD 177409 (HIP 94154), HD 144660 (HIP 
78983), HD 151528 (HIP 82260), HD 95338 (HIP 53719), HD 127352 (HIP 70973), HD 183063 
(HIP 95722), HD 105328 (HIP 59143), HD 204385 (HIP 106213), HD 145809 (HIP 79524), HD 
172675 (HIP 91645), HD 149013 (HIP 81010), HD 185615 (HIP 96854) and HD 138885 (HIP 
76292). 

The seven PlanetPol stars that define the upper end of the region RA = 16 H — 20 H polarization 
envelope ( ^5.31 ) are listed in Table [3] 
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C. Polarization vrs Distance relation 



Polarization of optical starlight in the interstellar medium is due to magnetically aligned 
dichroic dust grains, so that there is a correlation bet ween the upper envelope of polarization and 
interstellar extinc tion for a given set of observations (|Serkowski et al.lll975t Fosalba et al.l 12002 : 
Anders son! |20 1 2[) . The usual term "upper envelope" describes the relation between the maximum 
observed polarizations for a set of stars as a function of extinction, i.e. for a given extinction no 
polarization is found above the upper envelope value by definition. Figure [7] shows that an upper 
envelope is also present in the polarization-distance plot. The scatter of polarization strengths for 
a given extinction is due to the patchy nature of the ISM and magnetic turbulence. Likewise, the 
presence of an upper envelope to t he polarization-dis tance plot in Figure [TJlikely results from either 
the patchy nature of nearby ISM (IFrisch et al.ll201l[) . or magnetic turbulence, or both. 



The polarization data represented in Figure |3]were collected over a timespan of a half of a 
century, using a diverse set of instruments with unknown systematics. The reasonableness of the 
statement that polarization increases with distance has been tested for the stars shown in the figure 
using the nonparametric Kendall's tau measure of bivariate correlation. For the set of stars where 
Pol/dPol> 3, the rank correlation coefficiet is r = 0.86, where a value of one indicates that there 
is a perfect association of increasing polarization with increasing distance. Selecting only the most 
accurate data points, Pol/dPol> 5, yields r = 0.89. However for the nearest stars used in this 
analysis, < 40 pc, this association is not demonstrated (r = 0.14 P/dP> 3), an effect that we 
attribute both to influence of systematic uncertainties on weak polarization measurements and to 
the patchy nature of nearby ISM. In §5.3.11 we discuss a special nearby region where increasing 
interstellar polarization is shown to be associated with increasing distance for stars within 40 pc. 



D. Asymmetry of very local interstellar radiation field 



The asymmetry of the local interstellar radiation field is related to the overall opacity of the 
Local Bubble. The brightes t of the far-UV sources are located in the third and fourth galactic 
quadrants (£= 180° — 360°. krondhalekar et al. 1980 ) because of the distribution of massive stars 
with respect to the Local Bubble void. Those same stars are bright at the 13.6 eV ionization 
edge of H° so th at in the abs ence of interstellar opacity the t— 180° — 360° interval will be more 
highly ionized (|Frischll2010|) . This distribution of radiation sources explains the relation between 
polarization and the radiation field at each star that is displayed in Figure [9j 

Figure [12] shows the different levels of far UV radiation at 975 A between galactic quadrants 
I and H (£= 0° - 180°) and III and IV (£= 180° - 360°). The sum of the stellar fluxes are shown 
at spherical surfaces around the Sun with radii of 25 parsecs and 45 parsecs, and displayed using 
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an aitoff projection. The mapped fluxes represen t the sum of stellar fluxes at each position, based 
on the twenty-five brightest OB stars at 975 A (|Opal & Wellerlll984|) . Flux units are photons cm 2 
s 1 A -1 . The fluxes are plotted with the same color-coding in both figures. The minimum and 



-i 



maximum flux values on the color bar correspond to 4.0 x 10 4 and 8.0 x 10 4 cm" 2 s _1 A 
respectively. Low radiation fluxes in the region of large PlanetPol polarizations are the result of 
the locations of hot stars. 
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Table 1 . Various Directions 







Motive 




r] p o 




Direction of best-fitting Magnetic Field 






Polarization, Paper I, unweighted: 






Ecliptic 263 


37 


uncertainties ±35 


Galactic 37 


23 


uncertainties ±35 


Polarization, this Paper, unweighted: 






Ecliptic 263±\° 


37 ± 15 




Galactic 37 ± 15 


22 ± 15 




Polarization, this Paper, weighted: 






Ecliptic 263±^ 


47 ± 15 




Galactic 47 ± 20 


25 ±20 




ISMF from Center of Ribbon arc: 




[2] 


Ecliptic 221 ± 4 


39 ±4 


(see text) 


Galactic 33 ± 4 


55 ±4 




Upwind direction of interstellar He° flow through heliosphere 


V=23.2 ± 0.3 km s -1 , T= 6300 ± 390 K [1] 


Ecliptic 259.00 ±0.47 


4.98 ±0.21 




Galactic 5.25 ± 0.24 


12.03 ±0.51 




Quadrant III pulsars 






Ecliptic 232 


18 


[6] 


Galactic 5 


42 


[6] 


Heliotail in globally distributed IBEX ENA flux 




Ecliptic 30 ± 30 


0±30 


[5] 


Galactic 146 ± 30 


-49 ± 30 


[5] 


Direction of tail-in cosmic ray asymmetries 






Sub-TeV anisotropics: 






Ecliptic 90 


-47 


cone half-width=68 [3] 


Galactic 230 


-21 


cone half-width=68 [3] 


Ecliptic 66 


-36 


center of Gaussian fit [4] 


Galactic 211 


-35 


center of Gaussian fit [4] 



References. — HI iMcComas et alJ <l2012l) [2] iFunsten et alJ (l2009h [3] iNagashima et all (119981) [4] lHall et all (11999b [5] 
Schwadron etail d201 lb T61 ISalvatil (bold) 
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Table 2. New LNA and NOT Polarization Data 



HD 


I 


b 


Dist. 


Pol 


PAcol 


PAgal 


Source 




deg 


deg 


pc 


10~ 5 


deg 


deg 






1581 


309 


-53 


8 


15 ± 16 


73 ±52 


60 ±52 


LNA 


2025 


32 


-83 


18 


142 ±88 


125 ± 17 


34 ± 17 


LNA 


3443 


69 


-84 


15 


19 ± 13 


138 ±52 


2 ±52 


LNA 


4628 


120 


-57 


7 


7± 11 


86 ±52 


88 ±52 


LNA 


5133 


297 


-87 


14 


45 ±73 


75 ±52 


72 ±52 


LNA 


7570 


290 


-71 


15 


41 ±14 


70 ±9 


86 ±9 


LNA 


10360 


289 


-59 


8 


22 ± 11 


158 ± 14 


179 ± 14 


LNA 


14412 


209 


-70 


12 


8± 17 


28 ±52 


131 ±52 


LNA 


16160 


162 


-48 


7 


20 ± 15 


131 ±52 


95 ±52 


LNA 


23356 


210 


-50 


14 


449 ± 53 


7±3 


117 ± 3 


LNA 


125072 


313 


2 


11 


29 ± 135 


156 ±52 


174 ±52 


LNA 


129502 


346 


47 


18 


3±3 


110 ±24 


148 ± 24 


NOT 


130819 


340 


38 


23 


6±3 


81 ±14 


114± 14 


NOT 


131923 


323 


10 


24 


47 ±24 


41 ±14 


68 ± 14 


LNA 


131976 


338 


33 


7 


63 ±65 


30 ±52 


63 ±52 


LNA 


131977 


338 


33 


5 


55 ±24 


66 ± 12 


99 ± 12 


LNA 


134987 


339 


27 


25 


14 ±2 


65 ±4 


100 ±4 


NOT 


136894 


340 


24 


28 


14 ±4 


70 ±8 


106 ±8 


NOT 


141272 


369 


40 


21 


63 ± 32 


56 ± 13 


110 ± 13 


KVA 


144253 


352 


23 


18 


13 ±3 


69 ±7 


115 ±7 


NOT 


144585 


358 


27 


28 


2±2 


56 ±23 


104 ±23 


NOT 


152311 





14 


28 


3±3 


39 ±20 


91 ±20 


NOT 


153631 


7 


17 


26 


5±3 


6± 13 


61 ±13 


NOT 


154088 


355 


7 


18 


6±3 


60 ± 11 


113 ± 11 


NOT 


156897 


3 


8 


17 


5±3 


127 ± 13 


3± 13 


NOT 


160346 


27 


18 


10 


30 ±28 


173 ±52 


55 ±52 


LNA 


161096 


29 


17 


25 


5±0 


166 ±0 


49 ±0 


Wik 


165222 


24 


9 


7 


17 ± 18 


23 ±52 


84 ±52 


LNA 


167665 


4 


-5 


29 


4±3 


113 ± 16 


174 ± 16 


NOT 


169916 


7 


-6 


23 


7±3 


109 ± 12 


171 ±12 


NOT 


169916 


7 


-6 


24 


36 ± 7 


109 ±7 


172 ±7 


KVA 


170657 


13 


-4 


13 


4±3 


66 ± 18 


128 ± 18 


NOT 


172051 


12 


-6 


13 


2±3 


110 ±26 


173 ±26 


NOT 


176029 


38 


1 


11 


25 ±23 


141 ±52 


23 ±52 


LNA 


177716 


9 


-15 


37 


38 ± 7 


113 ± 7 


179 ±7 


KVA 
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Table 2 — Continued 



HD 


t 


b 


Dist. 


Pol 


PAcd 


PAgal 


Source 




deg 


deg 


pc 


10~ 5 


deg 


deg 




178428 


49 


4 


21 


2±3 


63 ±23 


125 ±23 


NOT 


179949 


13 


-15 


27 


10 ± 3 


77 ± 7 


144 ± 7 


NOT 


180409 


25 


-10 


28 


4± 2 


95 ± 16 


159 ± 16 


NOT 


180409 


27 


-9 


28 


23 ± 14 


50 ± 16 


113 ± 16 


LNA 


180617 


40 


-2 


5 


37 ±82 


25 ± 52 


87 ± 52 


LNA 


184489 


41 


-7 


14 


7± 16 


73 ±52 


135 ± 52 


LNA 


184509 


18 


-18 


31 


7± 3 


108 ± 11 


175 ± 11 


NOT 


184985 


25 


-16 


30 


8±3 


88 ±9 


153 ± 9 


NOT 


186427 


83 


13 


21 


5 ± 2 


111 ± 9 


174 ± 9 


NOT 


188088 


18 


-22 


14 


10 ± 12 


138 ±52 


26 ±52 


LNA 


189245 


8 


-27 


20 


38 ± 12 


144 ±9 


37 ±9 


LNA 


189340 


33 


-18 


24 


28 ±21 


166 ±52 


50 ±52 


LNA 


189931 


4 


-29 


27 


17 ± 16 


108 ±52 


3 ±52 


LNA 


191408 


5 


-30 


6 


15 ±35 


65 ±52 


141 ± 52 


LNA 


191849 


355 


-32 


6 


16 ±5 


167 ±8 


69 ±8 


LNA 


191862 


30 


-23 


28 


1±3 


112 ±34 


178 ± 34 


NOT 


192310 


15 


-28 


8 


36 ±21 


156 ±16 


48 ± 16 


LNA 


194640 


13 


-32 


19 


47 ± 1760 


17 ±52 


91 ± 52 


LNA 


194943 


26 


-29 


30 


1±3 


170 ± 35 


58 ±35 


NOT 


196761 


22 


-32 


14 


25 ± 18 


128 ±52 


19 ±52 


LNA 


202560 


6 


-43 


3 


10 ±32 


89 ±52 


175 ±52 


LNA 


207129 


351 


-48 


15 


32 ±9 


93 ±8 


14 ±8 


LNA 


211415 


340 


-51 


13 


13 ±5 


41 ± 10 


156 ± 10 


LNA 


213042 


18 


-58 


15 


18 ± 18 


108 ±52 


14 ±52 


LNA 


222237 


311 


-44 


11 


86 ±97 


52 ±52 


28 ±52 


LNA 



Table 3. PlanetPol* stars in polarization envelope 



HD 


hame 


Coordinates 


Distance 


Polarization 


PAra 






deg 


pc 


ppm 


deg 


172167 


a Lyr 


67.4, 19.2 


7.8 


17.2 ± 1.1 


34.5 ±1.4 


159561 


a Oph 


35.9, 22.6 


14.3 


23.4 ± 2.0 


30.8 ±2.4 


161868 


7 Oph 


28.0, 15.4 


29.1 


40.8 ±3.1 


28.5 ±2.1 


164058 


7 Dra 


79.1,29.2 


45.2 


73.3 ±1.2 


145.0 ±0.5 


186882 




78.7, 10.2 


52.4 


108.0 ±2.0 


22.8 ±0.5 


168775 




63.5 ,21.5 


70.4 


106.5 ±2.7 


6.1 ±0.7 


189319 




58.0 , -5.2 


84.0 


199.5 ± 1.4 


56.1 ±0.2 



Baile y et al. 



(120101) 
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Fig. 1. — The configuration of the ISM associated with Loop I is displayed. Filled contours show 
interstellar dust 100 within parsecs and the short lines show interstellar polarization directions. The 
extended gray -blue curved lines show the best-fitting ISMF for the region within 40 pc and 90° of 
the heliosphere nose ( §4.2[ Table [T]). Star distance is coded as follows: within 40 pc (cyan blue), 
40-100 pc (orange) and 100-150 pc (purple). Line length is unrelated to polarizat ion strength 



The Loop I ISMF direction is clearly defined by the polarizations of distant stars (e.g. ISantos et al. 



201 lh . but also appears to extend close to the Sun ( §5.31) . The lightest contours show maximum 



dust color excess values, ^(B-V), while the darkest regions with E(B-Y)< 0.01 mag show that 
the interior void of the Loop I bubble extends to the Sun location. The color excess E(B-V) 
contours are based on t he photometric and astrometric data for stars brigh ter than V =9 mag in 
the Hipparcos catalog (|Perrymanlll997f) and the absolute stellar colors in ICoxl (|2000|) : E(B-V) 
values are then smoothed spatially over ~ ±13° intervals for overlapping distances (considering 
uncertainties). Variable stars with 5V > 0.06 mag, as indicated by the Hipparcos variability index 
H6, have been omitted. The ^(B-V) contour levels of 0.01, 0.04, 0.09, 0.17 mag correspond to 
column densities \og(N(H°+H 2 )) of 19. 76, 20.37, 20.72, and 20.99 cm" 2 , for N(H°+H 2 )/E(B- 
V)= 5.8 x 10 21 atoms cm -2 mag— 1 (|Bohlin et al.1 1 1978b . Note the third and fourth quadrant 
voids (dark region between £~210°-0°) compared to the nearby brightest interstellar radiation field 
(Appendix |D]). The figure is centered on the galactic center, wi th longitude in creasing towards the 
left. The polarization data are from the sources listed in 03land lHeilesl (1200 1|) . 
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Fig. 2. — Column densities are plotted as a function of distance toward nearby stars within 40 pc 



and 9 0° of the heliosphere nose. The H° column density is found fr om either Lva data ([Wood et al 
20051) . or the the sum of D° components using D°/H°= 1.5 x 10" 5 dRedfield & Linskvll2004al) . The 



uncertainties are plotted as 20% of the column density. 
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Fig. 3. — Interstellar polarizations with P/dP> 3.0, and for stars within 90° of the heliosphere 
nose, are plotted as a function of distance. Data sources are color-coded, according to the scheme 
in the figure. Dar k purple points are the PlanetPol data. The boxes indicate stars in region RA = 
16 H — 20 H where bailey et al.l (|2010l) found that polarizations increase with distance. This region 
overlaps the North Polar Spur region, which forms the tangential sections of Loop I. PlanetPol 
requires bright stars to achieve high accuracy, which suggests that systematically lower PlanetPol 
polarizations beyond 40 pc (and perhaps closer) may be due to selection effects related to th e 
patchiness of foreground dust. Data are from Table|2] the references listed in ^3} and lHeilesI (|2000|) . 
These data were collected over a variety of wavelength ranges (©, and dif ferences in polarizatio n 
strengths will occur because of the wavelength dependence of polarization (|Serkowski et al. I fl975l) . 
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Fig. 4. — The polarization vectors of stars within ~ 40 pc are shown in both ecliptic and galactic 
coordinate systems, and color-coded for the data source, as indicat ed in the figure legend. The 



polarization data are from this paper (LNA and NOT, PlanetPol (Bai ley et al. 



I (LICK, KVA), ISantos et alJEoiI CTIO), and 20th century data dPiirolalll977 



201 0|, PP), Paper 



Tinbergenlfl982[ 



PT). Both plots are centered close to the heliosphere nose located at ecliptic coordinates (purple 
circle) of A=259°, (3=5°; longitude increases towards the left. Data with P/dP< 2.5 are plotted 
as dots, and vector lengths do not indicate polarization strength. The Comp ton- Getting corrected 
ENA fluxes at 1 keV are plotted for directions where the ENA count rates are larger than 113 
counts cm~ 2 s _1 sr _1 keV" 1 , which is y 1.5 times the mean ENA flux at 1 keV as measured by the 
IBEX-HI instrument ( McComas et al. 2009h. The grid of dotted lines show the ISMF determined 
from the center of the Ribbon arc (|Funsten et al. 2009 , Table 0). 



-44- 




Fig. 5. — The value of the function Fj = sm(6 ) n i ) for the unweighted fit (G n =l, eqn. [D is 
plotted over a regular grid of i possible interstellar field directions, based on unweighted position 
angles ( §4.11) . The function is color-coded and plotted in the ecliptic coordinates (left) and galactic 
coordinates (right), with the left plot centered on the ecliptic nose at A ~ 259°, and the right 
plot centered on the galactic center (16° from the ecliptic nose). The gray dotted grid shows the 
best-fitting ISMF direction for the unweighted fit listed at the top of the figures (Tabled). Stellar 
polarization vectors are shown superimposed on black dots for visual clarity. 
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Fig. 6. — The value of the function F[ = sin(9 nji ) /G n (eqn. Q]) evaluated over a regular grid 
of i possible interstellar field directions, based on weighted position angles with G n given by eqn. 
[2l The symbol size is coded to increase with polarization strength, and the dots represent stars 
where the observed polarization is not statistically significant. The function is color-coded and 
plotted in ecliptic (left) and galactic (right) coordinates. The gray dotted grid shows the best-fit 
ISMF direction given at the top of the figures (Tabled])- The secondary weaker minimum centered 
near the heliosphere nose, at plot centers, is dominated by the randomly distributed position angles 
contributed by low significance data points. See Figure |5]for additional information. 
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Fig. 7. — Right: The polarizations of the PlanetPol stars with P/dP> 3.0, distance D < 94 
pc, and inside the RA = 16 H — 20 H region are plotted as a function of distance (see text). The 
upper envelope of the increase of polarization with distance is shown by the fil led circles. This 



upper envelope is analogous to the increase of polarization with color excess (|Serkowski et al. 



1975b . A linear fit to the seven stars in the upper envelope (listed in Appendix |B]) yields logP 



1.200 + 0.013 * D, where P is in ppm and D is in parsecs. Left: The positions of these stars are 
plotted in galactic coordinates. 
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Fig. 8. — Left: Polarization position angles (PA R a) in equatorial coordinates are plotted against 
the star distance for the PlanetPol data in the RA = 16 H — 20 H region. Filled circles indicate the 
seven stars that form the upper envelope of the polarization envelope (Figure [7] right). Except for 
the outlier data point HD 164058 (at 45 pc, PA R a=145°, §5.31) . the polarization position angles of 
the envelope stars vary smoothly for the stars within 55 pc. A linear fit was made to the position 
angles of envelope stars within 55 pc, excluding the outlier HD 164058. The result is PAra= 
A - B • D pc , where A=35.97 ± 1.41, B= 0.252 ± 0.030, \ 2 = 0.559 and D pc is the distance in 
parsecs. The probability that the x 2 value is larger than 0.559 by chance is 0.76, indicating that a 
line provides a satisfactory fit to the position angles of the nearest envelope stars. A measure of 
magnetic turbulence is provided by the 23° standard deviation of the position angles of all stars 
within 55 pc, excluding the outlier, with respect to the linear fit to the nearest envelope stars. 
Right: Polarization position angles rotated into galactic coordinates are plotted at the locations of 
the envelope stars within 70 pc (except for the outlier). Data are labeled by the star distance (Table 
[3]). The best-fitting ISMF from the weighted fit (red box) and the center of the IBEX Ribbon arc 
(blue box) are plotted for comparison. Uncertainties on the position angles are negligible. Note 
that the position angle of the closest star HD 172167, 8 pc, is directed towards the best-fitting 
ISMF. 
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Fig. 9. — Stellar polarizations with 3a detections are plotted agai nst the FUV radiation field at each 
star for the 975 A fluxes from the 25 brightest stars measured by lOpal & Welled (|1 9 84I) . Dots show 
stars in the RA = 16 H — 20 H region, and diamonds show stars outside of this region. Black/gray 
symbols indicate stars within/beyond 45 pc. Symbol sizes are inversely proportional to the distance 
of the star. The seven envelope stars in Figure [7] are circled. The dashed/dotted lines show fits 
between the polarizations and 975 A fluxes for stars inside/outside of region RA = 16 H — 20 H . 
Polarization increases strongly at low FUV fluxes, but the effect is probably due to the relative 
locations of the most polarized distant stars in the first galactic quadrant, and the bright FUV 
sources that are mainly located in the third galactic quadrant (see Appendix ID)) . 
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Fig. 10. — Polarization strengths tend to increase with the ENA fluxes for the polarizations with the 
highest statistical significance. Left: The abscissa shows the polarization for stars within 40 pc and 
90° of the heliosphere nose, for P/a > 2.5 where P, a are the polarization and la uncertainty. The 
ordinate shows the mean IBEX-HI 1.1 keV ENA fluxes towards the same stars, where the IBEX 
fluxes are calculated as the average value for pixels that have S/N> 3.5 and are within 6 degrees 
of the star. The polarization sources are identified by color, with PT (black) representing 20th 
century data, PP (purple) representing PlanetPol data, and LNA, NOT, and LICK data representing 
data acquired for this study (blue, red, cyan, see The data shown in the left figure are divided 
into four groups, denoted by A, B, C, and D, given by (A) 1 x 1(T 5 < P < 11 x 10~ 5 , all 
ENA flux levels;(B) P < 1 x 1CT 5 , all ENA flux levels; (C) 11 x 1CT 5 < P < 40 x 1CT 5 and 
ENA fluxes > 130; (D) 11 x 10~ 5 < P < 40 x 10~ 5 and ENA fluxes < 130. The ENA fluxes 
increase with polarization strength for stars in group A. PlanetPol stars in group A are located in 
the RA = 16 H — 20 H region. Right: Stars in the four groups are overplotted on the IBEX Ribbon 
(1 keV), and color-coded according to the group in the left figure: red=A, B=blue, C=orange, and 
D=light-blue. Symbol size is proportional to polarization strength. The polarizations in groups A 
and B are recent data with lower uncertainties; group A in particular shows polarization increasing 
with ENA fluxes. The light-blue data points (group D) tend to be older data with large uncertainties. 
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Fig. 1 1 . — The great circle (gray line) in the sky that evenly separates the hot (red asterisk) and cold 
(blue asterisk) poles of the cosmic microwave background dipole anisotropy is plotted in ecliptic 
(left) and galactic (right coordinates). The best-fitting ISMF is shown in purple on each plot, 
with the dot/polygon giving the direction/uncertainties of the weighted fit in Table [Q The great 
circle equidistant between the CMB dipoles intersects both the best-fitting local ISMF direction 
and the heliosphere nose (red square, see Tabled), to within the uncertainties. The blue open circle 
gives the ISMF direction derived from the center of the IBEX Ribbon arc. The large blue circular 
region (light gray ) shows the location of the tail-in anisotropy for y 500 GCRs GeV mo deled by 



Hall et all (11999h . with the extent defined by the 68° width from iNagashima et all (119981) . The 
regions outlined by thick black lines show the observed heliotail that is deflected ~ 44° to the west 
of the do wnwind gas direction (red triangle), and is prominent in maps of the globally distributed 
ENAs (|Schwadron et al. 2011 ). Each plot is centered on the longitude of 0° in the respective 
coordinate system, with longitude increasing towards the left. 
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Fig. 12. — Diffuse interstellar fluxes at 975 A are plotted for spheres of radius 25 pc (left) and 
45 pc (right) around the Sun. The figure is centered on the galactic center with galactic longitude 
increasing towards the left. These fluxes are calculated fro m the combined fluxe s of the 25 brightest 
stars in the sky at 975 A based on the measurement s of Opal & Wellei dl984 ). Most of that flux 



originates in the third and fourth galactic quadrants (|Gondhalekar et al 



1980), which leads to the 



strong radiation fluxes in quadrants III and IV (right figure, £> 180°) when compared to weak 
fluxes in quadrants I and II (£< 180°). This asymmetry in the diffuse interstellar radiation field is 
due to the distribution of massive stars around the Local Bubble void. The more distant stars in the 
RA = 16 H — 20 H interval are exposed to a lower ambient interstellar radiation field than in regions 
closer to the Sun. Color bars are the same for each figure, and are plotted with arbitrary units that 
are scaled between the minimum and maximum fluxes of 4.0 x 10 4 and 8.0 x 10 4 photons cm~ 2 
s _1 A -1 , respectively. 



